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We propose to use a nonlinear resonator for projective readout, classical memory, and feedback 
for a superconducting qubit. Keeping the classical controller at cryogenic temperatures sidesteps 
many of the inefhciencies inherent in two-way communication between temperature stages in typical 
systems with room temperature controllers, and avoids increasing the cryogenic heat load. This 
controller may find a broad range of uses in multi-qubit systems, but here we analyze two specific 
demonstrative cases in single qubit-control. In the first case, the nonlinear controller is used to 
initialize the qubit in a definite eigenstate. And in the second case, the qubit’s state is read into the 
controller’s classical memory and used to reinstate the measured state after the qubit has decayed. 

We analyze the properties of this system and we show simulations of the time evolution for the full 
system dynamics. 

PACS numbers: 03.67.Lx, 42.50.Dv, 85.25.-j, 42.50.Pq 


I. INTRODUCTION 

The field of circuit QED [Ull], which integrates super¬ 
conducting resonators and qubits, has recently become a 
promising candidate for quantum information processing 
[3H5] due to dramatic improvements in the quantum tech¬ 
nologies [SHE]. Coherent operations can be performed fast 
and with high fidelities HIE], and different strategies for 
state initialization [TOUIS] and readout [HHn] have been 
demonstrated. However, the pursuit of large-scale quan¬ 
tum information processing is necessarily accompanied 
by a massive, further increase in quantum and classical 
system complexity for initialization, measurement, and 
control. As complexity increases, the inefficiencies inher¬ 
ent in the standard approach to measurement and control 
using external, room temperature apparatuses will be¬ 
come increasingly burdensome. Currently, circuit QED 
components are operated below 100 mK, and measure¬ 
ment signals used for readout and conditional control are 
sent through an expensive and meter-scale amplification 
chain up to room temperature. Similarly, controlling the 
quantum components requires transmitting signals from 
room temperature through ^60 dB of cold attenuation 
to the <100 mK-stage [5]. Thus, the current approach 
to inter-temperature, two-way communication between 
quantum and classical hardware involves inefhciencies in 
size, power, cost and speed that may prove incompat¬ 
ible with quantum information processing beyond the 
few-qubit regime. 

Another possible approach that sidesteps these ineffi¬ 
ciencies would be to perform much of the measurement 
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and actuation at the same temperature and with the same 
technology as the quantum hardware. In this article we 
analyze a specihc case of this alternate approach in which 
a nonlinear, superconducting microwave resonator |18j 
effectively reads-out and feeds back to a superconducting 
3d-qubit [^. The nonlinear resonator acts like a Joseph- 
son bifurcation amplffier (JBA) |15[ fl^ [20] in that the 
circuit’s Kerr-nonlinearity is driven and latches into one 
of two bistable states [H] , depending on the state of the 
qubit. In this way, the JBA implements a quantum non 
demolition (QND) readout of the qubit state and stores 
this state in its one-bit classical memory. Moreover, this 
classical memory-state may be maintained even if the 
qubit subsequently decays. However, the Kerr-resonator 
does not play the standard JBA role, as communicat¬ 
ing this state to the room temperature equipment (and 
experimenters) is not critical. Rather, a state stored 
in the Kerr-resonator is fed back directly to implement 
conditional qubit operations. 

Although we anticipate that natural advantages may be 
most practical in multi-qubit systems, here we analyze two, 
specific demonstrative cases in single qubit-control. In the 
first case, the nonlinear controller is used to initialize the 
qubit in a definite eigenstate. And in the second case, the 
qubit’s state is read into the controller’s classical memory, 
where it is stored for an indefinite period of time, and 
then used to reinstate the measured state after the qubit 
has decayed. The outline of of the article is as follows: 
In Sec. H, we introduce the system components and 
formalise the protocols for the two schemes. In Sec. HI, 
we recall the basic properties of the Josephson bifurcation 
amplifier, which is a central component in the proposed 
microwave network and we also present our theoretical 
formalism for the fully coupled system. In Sec. IV, 
we present numerical simulations with experimentally 
realistic component parameters and in Sec. V we compare 
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FIG. 1. (Color online) A qubit-cavity system is connected to 
a Kerr-resonator through a tee-junction. The qubit with a 
transition frequency Qqb is embedded in the 3d-cavity with 
a resonance at utb and the coupling port for the 3d-cavity 
is placed a distance dcav from the tee-junction. The Kerr- 
resonator, which consists of an array of SQUIDs that allows 
for a tunable frequency cOa, is placed a distance dx from the 
tee-junction. The Kerr-resonator is driven by a strong drive ad, 
while both the Kerr-resonator and the 3d-cavity are driven by 
the field Oin. Both drives are operated at the same frequency 


our scheme with other approaches for feedback and control. 
Sec. VI concludes the article and presents an outlook. 


off-resonant microwave signal, microwave power builds 
up in the mode and the signal is reflected with a certain 
phase shift. 

While the center frequency of the 3d-cavity is fixed, 
the Kerr-resonator’s center frequency is highly dynamic. 
Firstly, applying an external flux bias through the Kerr- 
resonator’s dc superconducting quantum interference de¬ 
vices (SQUIDs) via the flux control line will shift the 
center freauencv|18j. Secondly, the amount of power 
circulating inside the Kerr-resonator also decreases the 
center frequency. This power-induced shift is a highly non- 


(a) 

Simplified schematic 


3d cavity+qubit 


lb> 




Kerr resonator 


P 


(c) 

Kerr-resonator 

bifurcation 


(b) 

3d-cavity response 
On-resonance drive: 


Off-resonance drive; 


□ 


or 


(d) 

3d cavity and qubit: 


| 0 > 




\ 0^2 0^2“^ 

iai+a2p<|a,J^ |ai+a2r>|a,±p 


| 1 > 


© 


II. QUALITATIVE DESCRIPTION OF THE 
SCHEME 

The proposed setup is sketched in Fig. The two 
critical components are a two-port Kerr-resonator m and 
a 3d-cavity containing a superconducting qubit . These 
are interconnected via two transmission lines and a tee- 
junction. The microwave network has three connections 
to the external environment. The two most important 
are the transmission lines emanating from the bottom, 
weakly coupled port of the Kerr-resonator, and from the 
tee-port. A less critical port from the 3d-cavity uses a 
Purcell filter [niiinMi to control the photon number 
in the 3d-cavity and helps optimize performance. We 
will only discuss this port briefly in this section as the 
essential scheme can be understood without it. Dynamical 
control of the network is achieved by two microwave drives, 
labelled by the incoming amplitudes ad and and a 
flux control line that controls the center frequency of the 
Kerr-resonator. 

The basic operations that the network exploits are 
depicted schematically in Fig.J^ When a resonant mode, 
like an empty 3d-cavity in Fig.^(b), is driven by a slightly 


FIG. 2. (Golor online) (a) A simplified schematic of Fig. 
introducing symbols used in panels (b)-(d) and in Fig. The 
rectangle with a rectangular cut-out represents the 3d-cavity, 
the circle the qubit (in state \tp)), and the rectangle with a 
triangle cut-out the Kerr-resonator. Object colors are coarse 
representations of the center frequency of each component, 
(b) An empty 3d-cavity is driven by a resonant microwave 
drive (left-facing blue arrow), microwave power is built up 
in the 3d-cavity (cut-out is filled in blue), and the drive is 
reflected (right-facing arrow). When the same cavity is driven 
by a slightly-off- (light blue arrows) or far-off-resonant drive 
(green arrows), less or no power is built up inside the cavity. 
The reflected signal will now have a phase that depends on 
the sign of the detuned incident field (indicated by an either 
dotted or dashed left-facing arrow), (c) When the total power 
incident on the Kerr-resonator is sufficiently low the power 
built up is relatively low (lightly-filled triangle) but increases 
disproportionally when a certain pair of critical powers, |qc±P, 
are exceeded as detailed in the text (fully-filled triangle), (d) 
The qubit inside the 3d-cavity will shift the center frequency of 
the cavity by a small, but significant amount, as described in 
the text and yields a qubit-state dependent phase shift of the 
reflected signal. When the power is built up in the 3d-cavity, 
the center frequency of the qubit shifts proportionally (yellow 
circle accompanying the filled-in cut-out). 
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(a) JBA-like drive (b) Measurement 



(c) Stabilize memory 




(d) Feedback 




FIG. 3. (Color online) Protocol for measurement, classical memory storage, and feedback to a qubit, a) With the qubit in an 
initial state |i/)), the Kerr-resonator and drive are tuned to be co-resonant with the 3d-cavity (all three are colored blue), (b) 
Measurement operation as described in the text. (Note: the signal emitted by the Kerr-resonator is typically phase shifted as well, 
which we disregard in this section for simplicity.) (c) Stabilizing the classical memory state by tuning the Kerr-resonator and 
drive, and increasing drive power, (d) Using the classical memory state to apply a conditional 7r-pulse to the qubit, stabilizing it 
in the 11) state. 


linear function of the externally applied, near-resonant 
microwave drives (see Sec. Ill A), and Fig. [^(c) depicts 
the critical function for the purpose of this scheme. If the 
total incident power of red-detuned microwave drives is 
below a pair of detuning-dependent critical drive powers, 
lofcitP) the frequency shift will be negligible and little 
power is built up in the Kerr-resonator. If the incident 
power exceeds these critical values, though, the center 
frequency decreases significantly and it causes a stronger 
internal power to build up inside the resonator. If the 
incident power lies between the two critical values, the 
Kerr-resonator’s center frequency and internal power are 
bistable and hysteric (again, see Sec. Ill A for details). 


In Fig. (d) we depict the inclusion of a qubit, which 
changes the behavior of the 3d-cavity. The qubit state 
shifts the center frequency of the 3d-cavity and, conse¬ 
quently, the microwave signal from the 3d-cavity has one 
phase shift when the qubit is in the |0) state and another 
if it is in the |I) state. Conversely, the center frequency of 
the qubit is shifted in proportion to the amount of power 
built up in the 3d cavity. The details will be considered 
in Sec. IIIII 


Figure [^depicts a protocol that uses these functions to 
measure the qubit, store the result in a classical memory, 
and apply conditional feedback to stabilize the |I) state. 
In Fig.|^(a), the qubit is initially in an arbitrary state j'i/'), 
the Kerr-resonator is tuned using the flux control line to be 
near-resonant with the 3d-cavity, and the Kerr-resonator 
is driven through its bottom, weakly-coupled port. Some 
amount of drive power passes through the Kerr-resonator 


and drives the 3d-cavity, Fig.[^(b). The signal reflected 
by the 3d-cavity is partially lost in the remaining tee- 
port (this loss is not depicted) and partially reflected 
back to the Kerr-resonator. The phase of this reflected 
signal depends on the qubit state and interferes with 
the Kerr-resonator drive. Specifically, a |1) (|0)) qubit 
state will reflect a signal that interferes constructively 
(destructively) with the Kerr-resonator drive and pushes 
the total incident power above (below) the Kerr-resonators 
critical values. A disproportionally high (low) internal 
power then builds up in the Kerr-resonator. Thus, the 
system measures the qubit state and stores the result in 
the internal power level of the Kerr resonator. A high 
or low power level represents a 1 or 0 classical “memory” 
state, respectively. 


Figure [4] accompanies Fig. and depicts the detunings 
and amplitudes of the external controls used in the proto¬ 
col. The measurement just described occurs over a time 
Tm- During this interval, the Kerr resonator/3d-cavity 
detuning = uja—ujb is small, while the drive/3d-cavity 
detuning A'^ = uj^ — uib is zero with uia , Wb and ujd as the 
frequency of the Kerr-resonator, the 3d-cavity and the 
drive respectively. Furthermore, we depict in Fig. a 
drive applied at the Purcell-protected port, which adds 
additional signal of the qubit state to the Kerr-resonator 
but does not change the qualitative behaviour depicted 
in Fig. 3. Also shown in Fig. |^is the drive held of the 
Kerr-resonator, ad, which is slowly increased to a value 
between the two critical values, ac, between which the 
Kerr-resonator is bistable (see sec. IIIAl. 
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This memory state can be protected against subsequent 
changes in the qubit state by increasing the drive power 
and detuning the drive frequency and Kerr-resonator cen¬ 
ter frequency in tandem, as indicated by the changed color 
in Fig. (c). Although the phase of the drive reflected 
by the 3d-cavity is now almost independent of the qubit 
state, the total power incident on the Kerr-resonator re¬ 
mains relatively high or low depending on the most recent 
memory result. The classical memory state has been 
effectively stored and may be retained indefinitely. The 
external controls applied during stabilization is depicted 
in Fig. 1^ First, during a stabilization time Tg, the de¬ 
tuning of the Kerr-resonator and the drive amplitude are 
increased, which further separates the two critical values. 
This ensures that neither fluctuations in the field nor 
subsequent qubit dynamics will induce changes in the 


Measurement Stabilization Feedback 




FIG. 4. (Color online) Illustration of the external controls 
used in the feedback scheme. The upper panel sketches the 
power of the input drives, while the lower panel illustrates 
the detuning between the Kerr-resonator frequency and the 
3d-cavity, A), = uia—oJb (solid blue line) and between the drive 
and the 3d-cavity A); = ujd — uJb (dashed red line). During 
the first time-step, Tm, we apply a measurement driven by ad 
(solid red line) and an drive at the 3d-cavity, Qp (dashed gray 
line), which is followed by a stabilization time Ts, where the 
drive and Kerr-resonator parameters are changed to further 
stabilize the Kerr-resonator states by placing the drive far 
away from the critical input powers indicated by Oci (thinly 
dotted light red line). Simultanously the Kerr-resonator and 
the input-drive are tuned away from the 3d-cavity by At. After 
a waiting time Tmait in which no parameters are changed, we 
tune the Kerr-resonator back into resonance with the 3d-cavity. 
During the second Ta a classical input drive, atn (dashed blue 
line) is applied to the upper port and will interfere with the 
field from the Kerr-resonator. Finally, the 7r-pulse (dotted 
orange line) of duration T-^ excites the qubit state depending 
on the Kerr-resonator field strength. An optional pi-pulse 
(dotted green line) at the end of T.u,ait restores the qubit to its 
measured state, as described in the final paragraph of Sec. [11] 


classical memory state. To further decouple the classical 
memory from the qubit, the Kerr-resonator and drive 
are detuned from the 3d-cavity by a large amount, 
simultaneously with the stabilization. At this point, the 
memory may be held for an extended time T^ait- 

Finally, Fig. [^ (d) depicts how the qubit may be deter¬ 
ministically prepared in the |1) state, regardless of the 
initial state l"^). The Kerr-resonator and drive are tuned 
back toward the 3d-cavity’s frequency, resuming their 
values from the measurement step (Fig. [^ (b) ). This is 
depicted in Fig. [^ in which the detunings are brought back 
to their original values. If the Kerr-resonator memory 
state is 0, the power incident on the 3d-cavity through 
the Kerr-resonator is weak and can be eliminated by su¬ 
perimposing it with another weak and phase shifted drive 
through the remaining tee-port, depicted in Fig.[^(d) and 
by the dashed blue line in Fig. [^ In this case, no power is 
built up in the 3d-cavity and the qubit center frequency 
remains un-shifted. If the Kerr-resonator’s memory state 
is I, however, the drive incident on the 3d-cavity is rela¬ 
tively large and some power will build up in the cavity 
(despite the phase-shifted drive through the remaining 
tee-port) and shift the qubit center frequency. We are 
therefore left with a qubit whose center frequency is either 
shifted or not, depending on the classical memory state. 
If a third microwave pulse tuned to apply a 7r-rotation at 
the original qubit frequency is now applied at the tee-port 
(Fig.J^ (d) and the dotted orange line of duration in 
Fig. ff, a |0) qubit will be rotated to |1), but a |1) qubit 
will be off-resonant and will remain in |1). Thus, a qubit 
in the |1) state is prepared deterministically. 

The protocol closes an autonomous feedback loop at 
cryogenic temperatures but from the point of view of an 
operator at room temperature, it constitutes an open 
loop protocol. This is conceptually similar to quantum 
reservoir engineering schemes [HIISHIT], in the way that 
it changes the effective ground state and does not rely 
upon any conditional room temperature operations com¬ 
mon to conventional feedback schemes O IHlEHl Eg. 
The protocol, however, relies upon the performance of 
a coupling to a dissipative meter system and, thus, it 
falls conceptually in-between conventional feedback and 
quantum reservoir engineering schemes. 

We can further demonstrate that the classical memory 
may be usefully retained for an arbitrary amount of time. 
For example, the setup proposed here can return the clas¬ 
sical state stored in the memory to the qubit, even after 
the qubit has undergone Ti decay to |0) (and assuming no 
spontaneous excitation to |1) for demonstration purposes), 
i.e. Tujait ^ Ti. Just before we bring the Kerr-resonator 
back onto resonance, we apply a 7r-pulse on the bare qubit 
frequency bringing the qubit into the excited state, dotted 
green line in Fig.|^ After this pulse, we apply the same 
input field and tuning sequence as above. A resonant 
TT-pulse will now prepare the qubit in the excited state if 
and only if that was the state measured earlier. 










































III. QUANTITATIVE ANALYSIS OF THE FULL 
SYSTEM 
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The circuit illustrated in Fig. consists of different 
components. We shall now provide the formalism for the 
quantitative description of the systems and their dynamics. 
We will begin with a breif review of the readout mechanism 
followed by a complete analysis of the full system. 


A. The Josephson bifurcation readout 


Since the readout is fundamental for the operation of 
the scheme we will briefly review how a Kerr-resonator 
performs a projective readout of a qubit. A Kerr-resonator 
may be implemented as an array of SQUIDS in paral¬ 
lel with a capacitance jl8]. The SQUIDS provide the 
inductance for the equivalent LC-circuit and they per¬ 
mit tuning of the resonance of the Kerr-resonator by a 
magnetic flux. While a single SQUID provides a Kerr- 
nonlinear behaviour, an array will dilute the non-linearity 
and allow a much higher photon number in the Kerr- 
resonator [3D]. The resonator loss rates associated with 
the two coupling capacitors in Fig. [^are Kd and Ka and 
we assume a classical input field ad coupled with strength 
to the resonator mode. In a frame rotating with the 
drive frequency, the classical field amplitude in the 
resonator solves the non-linear equation. 


= -*A„ 


-iK\a\^a-'^^^^a-^dad, (I) 


where K denotes the Kerr-nonlinearity, which originates 
from the non-linear inductance of the SQUID-array. 

In the classical field approximation, the driven Kerr- 
resonator relaxes into a steady-state which is found by 
setting d = 0 in Eq. Q. This yields an equation for ad 
in terms of a. We are, at present, not concerned with the 
phase of the solution, and by multiplying the expression 
for ad with its complex conjugate, we find the non-linear 
relation between the driving field and the intracavity field, 

Kd\ad\^ = -b 2AaK\af 

+ A^|ar+ + (2) 


This expression yields unique values of for any a, 
while multiple values of the intra-cavity field strength 
a may be compatible with the same driving held. We 
analyze this possibility by solving = 0, which yields 

the critical photon numbers nc± = |q!c±Pi 


__ 2 ^/ i{Al + {Ka+Kd?/^) \ 

3A2 4A2 )■ 

(3) 

If Aa = A(j c = \/3/4:{Ka + Kd), we obtain a point of 
inflection in the equation ([^, while if we increase A^ 



FIG. 5. (Color online) Solution of the classical field equation 
0 with drive amplitudes slightly above and slightly below 
the critical drive, Oci- The inset shows how the two applied 
drive amplitudes linearly approach the critical value at time 
20/{Ka + Kd) indicated by the vertical dotted (green) line. A 
constant detuning Aa = 1.75Aa,c is assumed between the 
drive and the Kerr-resonator and the non-linearity is set to 
K = —0.012Aa. For the slightly higher drive power, we obtain 
a steady state with ~110 photons, which is much larger than 
the ~35 photons obtained for the only slightly weaker inital 
drive. 


beyond Aa,c we obtain two unique solutions, nc±- For 
driving strengths between the two corresponding field 
values, Eq. 0) has one unstable and two stable solutions. 
The Kerr-resonator bifurcated in the same way as a JBA, 
which derives its name from this bifurcation, and now 
a very weak initial perturbation is sufficient to make 
the resonator field attain different stationary solutions 
that depend on this perturbation. The Kerr-resonator 
will latch into one of these two solutions even when the 
perturbation is switched off. In Fig. [^we see an example 
of how this bifurcation manifests itself as a very large 
difference in photon number with only a small difference 
in input power. This large difference persists even when 
the input power is levelled out. The perturbation could be 
introduced in the upper input port to the Kerr-resonator, 
for example, from the 3d-cavity controlled by the qubit 
state in Fig.[^ In this case, a high fidelity QND detection 
of the qubit state results if the output fields are amplified 
and transferred to room temperature electronics HU HD], 
In this work, however, rather than transferring the output 
to room temperature, we will use the fact that the steady 
states of the Kerr-resonator are already classical states 
with high fidelity correlation with the qubit state. 

The time scale for readout is naturally related to the 
time of bifurcation, which in general is complicated to 
predict but is associated with 1/K. If the Kerr-resonator 
can bifurcate on a very short time scale, the readout 
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fidelity might be poor. This is because the strong coupling 
between the Kerr-resonator and the 3d-cavity, needed for 
fast readout, will lead to a larger Purcell-decay of the 
qubit and a very high photon number inside the 3d-cavity, 
which will cause the measurement to be non-QND. This 
obstacle can be mitigated by the Purcell-filtered port of 
Fig.[TJ allowing the cavity field to decay but maintaining 
the qubit lifetime. We can, thus, achieve very strong 
coupling and still remain in the QND-regime by effectively 
increasing the cavity decay rate. Eventually, though, the 
readout-speed and fidelity will still be limited by the 
dispersive coupling between the 3d-cavity and the qubit. 


B. Circuit analysis 


The Kerr-resonator and the qubit in the 3d-cavity are 
described by the Hamiltonians [301 El] : 


K 

Hk = fiUJaO^a -I- H—a'aa'a, (4) 

Hcav = f>^bb'^b + h^^az+hxbHaz, (5) 


with resonator eigenfrequencies uJa,b, a qubit frequency 
Hqf, and a dispersive interaction between the qubit and 
the 3d-cavity field with the strength y. Here az denotes 
the qubit operator such that ctj |1) = |1) and crz|0) = —10). 
The operators a (a^) and b (b^) describes the annihilation 
(creation) operators for the Kerr-resonator and the 3d- 
cavity respectively. 

The Kerr-resonator and the qubit-cavity are connected 
by one-dimensional waveguides. The position dependent 
phase operator, (time integral of the voltage) and 

charge density, q{x,t), along the waveguides obey wave- 
equations and they can be expanded on left and right 
propagating wave solutions |3H |33] , 


-b fi.c. (6) 

V 47 r 7 o 

^ I h7 ^ 

~q{x,t) = —iy J dcj Vw -b/i.c. (7) 

such that 


(l){x,t) = (j^ixjt) + cj){x,t). (8) 

In the above equations we have assumed that the wave¬ 
guide impedance Zq is constant and the same for all lines. 
We choose the origin of the position variable to be at 
the tee-junction. The left (right) arrow denotes the field 
moving towards (away from) the 3d-cavity and down (up) 
towards (away from) the Kerr-resonator. 

The waveguides are joined by a tee-junction, and Kir- 
choff’s equations yield [M] 

(l)out{t) = ^ ^ t) (9) 


for the reflected component, (jiout, in terms of the waves 
incident on the junction, (pin- Now (j)out{t) depends, in 

particular, on the fields from the Kerr-resonator, and 

from the 3d-cavity, (pcav 

Let qh denote the charge degree of freedom at the 3d- 
cavity and similarly let qa denote the charge at the Kerr- 
resonator, such that we obtain the boundary conditions 

m 


4>KidK,t) = + ZoPaqait), (10) 

- y ^ 

^ cav(.^cav') 4^cav(,^cavi t) + Zo(3i,qb{t), (11) 

where dx and dcav denote the lengths of the waveguides 
from the junction to the Kerr-resonator and the 3D cavity 
respectively. The quantities Pa and Pb are unit-less geo¬ 
metric constants that depend on the relation between the 
effective coupling capacitance and the total capacitance of 
each subsystem. The charge variables qa and qb, are given 
by the cavity mode operators (a, a^) and {b,b^), respec¬ 
tively. Through the so-called black-box quantization |3T] 
we assume the Hamiltonian of the qubit-cavity system 
to be of the form Eq. ([^. Following this description, qb 
now also includes a contribution from the qubit and this 
exactly leads to the well-known Purcell-decay of the qubit. 
We will omit the qubit dynamics for now, but it will be 
considered in Sec. lYl 

An input-output relation similar to Eq. © can also be 
written for the lower output port of the 3d-cavity shown 
in Fig 12 However, unlike the upper port, we will attach 
a Purcell filter [TBl to the output port. The sole 

purpose of this additional decay-channel for the 3d-cavity 
is to allow a very strong interaction between the Kerr- 
resonator and the 3d-cavity without exciting the 3d-cavity 
beyond the dispersive regime |35j . The extra cavity decay 
does therefore not by itself improve the qubit readout, 
but it allows for a much stronger coupling between the 
Kerr-resonator and the 3d-cavity. So as long as the total 
decay rate of the 3d-cavity remains in the same order of 
magnitude as the dispersive coupling to the qubit, we can 
increase the coupling to the Kerr-resonator by increasing 
the decay through the Purcell-fitered channel by the same 
amount. In its simplest form the Purcell filter is made 
by a transmission line stub connected to the output line. 
This stub acts as a A/4 impedance transducer, but only 
at a single frequency, w/. Now, if w/ = flqb, the filter will 
short out the environment of the transmission line and 
thus create an effectively vanishing density of states in the 
environment at the qubit frequency, leading to vanishing 
Purcell-decay of the qubit through the filtered channel. 
The 3d-cavity, however, operates at a different frequency 
and is allowed to decay through the filtered channel. 

We will assume a narrow frequency band of the trav¬ 
elling fields, kai « k, and we restrict the analysis to the 
situation of wave propagation times much shorter than 
other dynamical time scales of the systems. This allows 
elimination of the wave-guide observables, 4>{x, t) and 
q{x,t), and we obtain the Heisenberg-Langevin equations 
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for the a and b operators 


/3, 


. . . f “I” 

a = —lUJad — iKa' aa - - - a 

+ - y/l^ad, (12) 

b = -iuJbh - ixcTzb - ^ 6 

+ - y/t^ap. (13) 


a = —iAaa — iK\a\'^a — ^ + y/KaKbl^e''^^'’ 

y/ (1^) 

/3 = -iA{,/3 - ixWz)P - _|_ 

^ KpQ-p + ^^(t), (1"^) 


where the 3d-cavity and the Kerr resonator are coupled 
with the strength y/KaKb- In Eqs. (12) and (13), we have 
also included an extra input and decay channel for the 
Kerr resonator corresponding to the lower input port of 
Fig.l and for the 3d-cavity through the upper port (nb) 
and the Purcell filtered lower port (Kp). The operators 
ain and ad are the input field operators, while the phases 
4>a = dxk and = dcavk come from the propagation 
through the transmission lines to the Kerr-resonator and 
the 3d-cavity. 


Having Eqs. (12) and (131 we can in principle solve 


the full dynamics of the system, but since a and b are 
operators and the mean value of a might be on the order 
of ~10, we would need to expand the Hilbert space on 
a large number of Fock states to solve the full quantum 
dynamics. We, therefore, employ a stochastic mean-field 
description by making the c-number substitutions: 


ae 


— iuidt 




— ilddt 


b- 

O-d 


-t (cTz). 


For convenience we have transformed our equations to a 
frame rotating at the frequency of the external drive, oJd, 
which will make the detunings = uja — ujd and Af, = 
ujb — OJd appear in the equations of motion. We can make 
the complex field variables represent random samplings 
of the quantum Wigner quasiprobability distributions if 
we include stochastic terms, representing the quantum 
noise associated with the dissipation terms in (12) and 

ISS]. 

Following the prescription in |36j , we obtain fluctuating 
complex noise contributions 


Z,ait) = + Kd)iN + ^)dWa,l 

+ ^d)(W -|- —)dWa^2, (14) 

ib{t) = ^^p)i^ 4" 2 '^dWb^i 

+ i\j 4" ^p)i^ + 2^dWb^2, (15) 

where dWi^n denote independent Wiener processes with 
(dWi^n) = 0 and {dWf^) = dt and N is the mean photon 
number in the thermal bath coupled to the resonators. 

Adding these terms to the deterministic mean value 
equations, we obtain the equations of motion for a and 


The replacement of the full quantum operator equations 
by noisy mean field equations has been succesfully applied 
in a wide range of quantum optics problems, and due to 
the weak Kerr non-linearity we expect this description to 
correctly describe the mean value of the resonator field as 
long as the qubit is in one of the eigenstates of tr^. We refer 
to [37] for a recent comparison of full quantum and noisy 
mean field calculations in an optical network system with 
Kerr-nonlinearities. This system also shows bistability 
and switching dynamics, similar to ours and shows a 
fully satisfactory agreement between the methods. In the 
appendix we verify this agreement for simple instances of 
our problems, with the only discrepancies concerning the 
quantitative dynamics just at the bifurcation threshold, 
a minor concern, since the precise properties have to be 
calibrated in the experiments. 

As one sees from Eq. (17), the qubit state influences 
the driving of the 3d-cavity, which in turn affects the 
driving of the Kerr-resonator. In the regime where the 
qubit-dependent feedback drives the Kerr-resonator into 
two different stable states, this will lead to a projective 
measurement of the qubit, and storage of the outcome 
in the Kerr-resonator. The tunability of both A^, otd 
and ain now enables the implementation of the schemes 
described in Sec. ttn 


IV. NUMERICAL SIMULATIONS 


So far we have treated the qubit as a semi-classical 
object with two states, |0) and |1), each of which shifts 
the frequency of the 3d-cavity. In reality the qubit is 
described by the quantum mechanical density matrix, 
Pqh = Si j=o 1 il® dynamics are described 

by a master equation. 

However, we choose to describe the Kerr-resonator as 
well as the 3d-cavity by stochastic complex fields, where a 
single trajectory is determined by the particular instance 
of dW(t) restricting ourselves to classical-like states. Sim¬ 
ilarly we can unravel the master equation for the qubit 
into a non-Hermitian Schrodinger equation [SS] [35] 


zL 

dt 




A 


-l- 


qb 


■2xl/3p 




(Hd(t)cr+ + HJ(t)cr-)|V'), 


(18) 


where Aqb = Hqb — oJd and fid{t) describe a drive at the 
qubit frequency. The non-Hermitian Schrodinger equation 
is designed such that the norm of the wavefunction, lip) 























decays, due to the excited state decay rate 7 = l/Ti. 
When the norm becomes smaller than a random number 
R, drawn uniformly between 0 and 1, we apply the jump 
operator (t_ to the state and then we renormalize l"^). 
This procedure reproduces the average density matrix 
evolution of the qubit. 

A superposition state of the qubit generally yields ex¬ 
pectation values (az) = {ip\az\'4’) that differ from the 
eigenvalues ± 1 , but after the bifurcation readout of 
the qubit, we expect the state of the coupled system 
to be in a completely mixed state with different field 
states pertaining to the excited and ground state qubit, 
p = Po|0)(0| (g) pap -b Pi|l)(l| g Pa,i, with Pa,0(1) being 
the 0 (1) memory state of the Kerr-resonator. Here P„ 
denote the probabilities for the qubit to be in state 0 and 
1, and due to the QND interaction, these probabilities do 
not change during the readout. Thus, we can simulate 
the two classical qubit cases separately, and on average 
generate the full master equation evolution. 


Now, we can investigate the performance of the pro¬ 
posed scheme by numerically solving Eqs. 0 

and (18), but with the qubit in an eigenstate as de¬ 
scribed above leaving (cr^) = ±1. For the numer¬ 
ical simulation we need realistic parameters, so we 
choose a detuning between the qubit and the res¬ 
onator at Aqb = 27r X 1.2 GHz. Furthermore we as¬ 
sume that the coupling between the two systems are 
at p = 27r X 122 MHz, which yields a dispersive shift 
at X = g‘^Ec/{fiAqi,{Aqb — Ec/h)) = 27r x —2.5 MHz [ID] 
with the anharmonicity of the qubit given by 
Ec/h = 27 r X 300 MHz. We set the decay-rate for the 
3d-cavity to Kf, = 27r x 1 MHz and Kp = 27r x 4 MHz. 
For the Kerr-resonator we assume AT = 27r x (—0.4) MHz, 
Ka = 27r X 5 MHz and = 2Tr x 0.3 MHz with a de¬ 
tuning of the center frequency from the drive initially at 
Aa = 3.5 X (kc + Kd)- Finally, for the qubit we assume 

the intrinsic decay-rate 7^6 = 27r x 5 kHz and from the 

2 

cavity we get an additional Purcell-rate of jp = 

^qb 

which yields 7 = 7 qt, -I- 7 p = 27r x 15.3 kHz corresponding 
to a qubit lifetime of Ti = 10 ps. 


For the bifurcation readout our goal is a large difference 
in the Kerr-resonator photon number conditioned upon 
the qubit state. We find that a measurement time of 
Tm = 400 ns is sufficient to obtain a significant difference 
between the signals for the different qubit states compared 
to the noise of the input field. A shorter measurement time 
would be dominated by noise and yield a poor readout 
fidelity. During the first 80 ns, we slowly increase the 
drive towards ad = 33.02ydcd, which is just below the 
larger ac± set by Eq. (f^ . The Purcell protected port was 
driven by ypE^ap / i/l'Ki) = 8 MHz, with the i indicating a 
7r/2- phase shift of the drive compared to the rest of the 
drives. The next step is the stabilization where we linearly 
increase Aa by a factor of 1.7 and ad by a factor of 1.8. 
This increases the difference in amplitude between the 
two solutions of Eq. (§ and puts ad further away from 
both - as a consequence the probability of a spontaneous 


jump between the two solutions vanishes. To ensure that 
we remain on the right branch of the bifurcation curve 
Ts cannot be too small and we choose Tg = 150 ns. The 
stabilization time Tg might be chosen smaller if we apply 
a non-linear tuning scheme for Aa and ad, but with the 
linear increase we numerically find that smaller times 
introduce significant errors. 

The next step is the frequency tuning of the Kerr- 
resonator along with the drive. To ensure that we are 
indeed very far-detuned from the 3d-cavity after the de¬ 
tuning, we choose At = 27r x 30 MHz and to avoid exciting 
the Kerr-resonator while tuning the frequency, we cannot 
exceed dAa/dt = with k = Ka + Kd- This limits the 
detuning time to Ta = 100 ns. This detuning between 
the Kerr-resonator and the 3d-cavity is also occurring 
during Tg simultaneous with the increase in Aa (see Fig. 
1^. The waiting time T^ait can be chosen according to 
our applications. We then wait for a driving time Td = 15 
ns before the 7 r-pulse with Rabi-frequency = 27r x 7 
MHz and = 7r/(2nd) = 35.7 ns. 


A. State preparation and stabilization 

For the state preparation scheme we want to have 
the waiting time, T^uait, as small as possible in order to 
quickly initialize the state into the excited state. Since 
l /(«^6 + Kp) = 31 ns, the 3d-cavity is (almost) empty after 
the two tuning steps, so we can choose a very small waiting 
time T^ait = 25 ns. After T^uait we tune the resonator 
and drive frequency back again and we also apply the 
input field a^, which is of equal amplitude but opposite 
phase of the 0 state of the Kerr-resonator obtained when 
the qubit was measured to be in its ground state. When 
we are back on resonance such that uid = ujb, we apply 
the TT-pulse and we obtain numerically a pulse fidelity of 
98%. The imperfection arises from the field fluctuations 
in the 3d-cavity and due to the fluctuating field emitted 
from the Kerr-resonator. These fluctuations are the same 
fluctuations that limit the readout hdelity. We apply 
square 7r-pulses in our simulations, and a study of pulse¬ 
shaping is beyond the scope of the present work, but 
we note that more elaborate pulse shapes may prevent 
leakage of population into higher excited states of the 
weakly anharmonic transmon qubits m and to mitigate 
some types of quantum fluctuations [42] . 

In Fig. [^we show the evolution of and from the 
numerical simulation we achieve |/3p « OAucrit during 
the readout, where Ucrit = A^/(4gi^) is the critical photon 
number at which the dispersive regime breaks down |35| . 
The effective detuning of the qubit due to higher order 
effects is reduced by around 10 % at OAucrit (included in 
numerical calculations). We cannot improve the scheme 
by increasing the drive power, as we are already close 
to the maximum photon number where the dispersive 
coupling offers a valid description. We even notice that 
just after the bifurcation, the photon number is momen¬ 
tarily increased beyond the dispersive regime. While this 
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FIG. 6 . (Color online) Numerical simulations of the state 
preparation with parameters described in the text without 
qubit decay. In (a) we start with the qubit in |0) and bring 
the state into |1) with 95.7% fidelity while in (b) we start 
with the qubit in the |l)-state, where it remains with 98.5% 
Hdelity. The upper panel of each hgure displays the mean 
value of CTz, the middle panel displays the photon population in 
the 3d-cavity and the lower panel the population in the Kerr- 
resonator. In each panel we plot both the mean values for 50 
trajectories (thin lines) and the average over 200 trajectories 
(thick line). The vertical dotted lines indicate the different 
time intervals similar to Fig. |^and explained in the text. 

is not important for the readout since we have already 
transferred the state to the resonator, we need to take into 
account the so-called dressed dephasing |33], which adds 
an additional decoherence channel for the qubit propor¬ 
tional to the photon number in the 3d-cavity. This effect 
will effectively decrease the lifetime of the qubit state 
and since the typical time where each trajectory remains 
above ricrit is 25 ns, we expect a l%-2% additional error 
from the additional dephasing. In [43] it is also shown 
that above a saturation point around O.Sricriti dispersive 
readout cannot be improved, however higher power leads 
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FIG. 7. (Golor online) Stabilization of the qubit by running 
the same scheme as in Fig. j^but 20 times after each other. 
The parameters are all listed in the text. The dashed (red) 
line is the average of 400 trajectories, while the dotted (light 
red) line is an exponential Purcell-enhanced decay with the 
decay-rate of the qubit, I /7 = 10 ps. The solid (blue) curve 
shows the same but with a decay-rate of I /7 = 30 ps acheived 
with an additional Purcell-filter as explained in the text and 
the densely dotted (light blue) curve shows the corresponding 
exponential decay. 


to a faster bifurcation and more stability of the bifurcated 
states, so the optimal power of our scheme is expected 
to be much higher than just the optimal power for pure 
dispersive readout. If there is still a demand for lower 
photon number, around 0.1-0.2ncrit, our scheme will still 
work, but with a hdelity of only around 85%-90%. 

Figure demonstrates the ability of our scheme to 
deterministically prepare the qubit in the exited state 
from any unknown state. For simplicity we have not 
included qubit decay in these simulations. Starting in the 
excited state yields a succes probability of 98.1% and if 
the qubit starts in the ground start we achieve a success 
probability of 96.1%. The error originates partly from 
the imperfect 7r-pulse and partly from the fact that the 
readout-hdelity is only around 98%. 

Qubit decay is included with a hnite qubit lifetimes of 
10 fis (red line) and 30 /xs (blue line) in the simulations 
shown for 20 runs in Fig. [^with the same parameters 
as in Fig. The time averaged state hdelity is 89.9% 
and it can be increased to 93.5% by averaging only at the 
instants of the protocol, where the preparation scheme is 
completed and the hdelity is known to be highest. The 
main limitation for the stabilization is the assumed qubit 
lifetime. For example, installing an additional Purcell- 
hlter on the upper output port (possibly between the 3d- 
cavity and the tee-junction) increases the qubit lifetime 
to around 30 ^s and leads to an average hdelity with the 
excited state of 95.4% (and 96.8% by only averaging after 
each protocol ends) in our calculations. 
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FIG. 8. (Color online) Demonstration of the memory recovery 
scheme with parameters described in the text. The main figure 
shows the mean value of the qubit operator Oz averaged over 
200 trajectories for each. The solid light blue line assumes |1) 
as the initial state, while the dashed red line starts in |0). The 
TT-pulse is applied after 78 /rs, where we obtain a fidelity with 
the original state of 98.6% for |1) and 96.0% for |0). The upper 
inset shows Kerr-resonator excitations for 20 randomly chosen 
realizations and we see two clearly distinct stable states. The 
lower inset displays a zoom of the region around the recovery 
of the state. 


B. Qubit readout and memory 


To demonstrate that the classical memory state can be 
stored and reused in a feedback scheme at a later point in 
time, we will bring back the measured qubit state after a 
very long time. This scheme relies upon qubit decay into 
the ground state, but in practice qubits within a 3d-cavity 
tend to have a residual thermal excited-state population 
beyond what is expected from the cryogenic temperature 
pi] . This analysis is thus an idealisation to demonstrate 
the memory properties of the Kerr-resonator. We now 
apply a very long T^ait, and just before applying the 
tuning sequence, we apply a rr-pulse to bring the qubit 
into the excited state. The final steps are then the same 
as in the state-preparation scheme. We have chosen Ty^ait, 
such that we start the tuning after 78 ^s. In Fig. [^we 
illustrate the performance of this memory scheme. From 
400 trajectories, we obtain a fidelity for reading out |1) 
and bringing the qubit back into |1) of 98.6%, consistent 
with the state preparation scheme for the excited state. 
Similarly for |0) we find a fidelity of 96.0%. The memory 
scheme is limited by the imperfect 7r-pulses as discussed 
in the previous subsection, thus we have a smaller fidelity 
for the |0) state. There is also a small (<1%) remnant 
excitation from the qubit after Ty,ait- 


V. COMPARISON WITH OTHER 
APPROACHES TO CLOSED-LOOP QUBIT 
CONTROL 


As shown by the numerical simulations, the setup in 
Fig. [l] implements a feedback loop, but in general there is 
not a single correct way to implement closed-loop control 
of superconducting microwave qubits, and all-cryogenic 
and inter-temperature closed loop feedback systems are 
likely to be both applied in a more mature technology. 
Thus, it is important to clarify when all-cryogenic con¬ 
trol might be naturally competitive. Room temperature 
measurement-based control [laillilHllig benefits from 
mature, flexible, and commercially-available hardware 
and is easier to debug. But, as already mentioned, room 
temperature-based feedback control requires interfacing 
very different and remote technologies. As such, amplifiers 
must be state of the art in terms of simultaneous low-noise, 
high bandwidth, and gain, and gross energy inefficiency is 
unavoidable between qubits operating at fW power scales 
and Watt-scale room temperature equipment. In contrast, 
the components of our proposed scheme are proximate in 
a cold, low-noise environment, so exemplary amplification 
is less critical. Moreover, the qubit and controller feature 
“only” a factor of ~ 10^ difference in operating power, and 
further reduction in operation power is conceivable |45j . 
A general concern with cryogenic classical control is the 
heat dissipation fundamental to digital computation, but 
in the scheme here, dissipation occurs off-chip, as remote 
from the sensitive qubit and Kerr-resonator as necessary 
(unlike, for example, rapid single flux quantum controllers 

gSlIlT]). 

There is also no single correct approach to all-cryogenic 
feedback control of superconducting qubits. Often, these 
schemes are lumped together under the category of 
“autonomous control”, but there are important distinc¬ 
tions within this category as well. The most popular 
autonomous approaches today in superconducting mi¬ 
crowave systems are continuous pumping schemes, in¬ 
spired by quantum optics techniques [laiMsz]. Here, 
microwave transitions in a single component (usually a 
circuit QED qubit or qubits) are designed so that a con¬ 
tinuous microwave pump drives the qubit (s) into a desired 
state. These schemes are often understood in terms of 
an internal feedback loop in which a cavity and the qubit 
states respond to each other. We distinguish this from 
our approach in that the Kerr-resonator controller may be 
separated from the qubit without a fundamental change to 
the qubit’s operation. In contrast, remove the cavity from 
a circuit QED qubit and it becomes a very different device. 
As a consequence, our qubit and Kerr-resonator controller 
may be optimized with greater independence. Also, as 
the Kerr-resonator’s role is well-compartmentalized to a) 
measuring a digital microwave signal and b) performing 
feedback with conditional filtering of a microwave pulse, 
it is easy to see how this concept may apply more gen¬ 
erally [48] in superconducting microwave systems (See 
also discussion in Sec. VII. In contrast, strongly coupling 
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a general quantum circuit to a cavity leads to a highly 
perturbed energy spectrum and it is much less clear how 
to design such systems for general continuous pumping 
schemes. 


VI. DISCUSSION AND OUTLOOK 

We have proposed a new type of feedback scheme that 
aims at closing the feedback loop within the cryogenic 
environment of circuit QED experiments. To achieve 
this we connect a Kerr-resonator to a 3d-cavity with a 
qubit. The goal is to control the qubit without having 
any knowledge of the qubit state outside the cryostat. 

As an application of this type of feedback we have 
simulated three applications. At first we showed that we 
can prepare the qubit in the excited state with an average 
fidelity of 97.1%. Including qubit decay we showed that we 
can continuously stabilize the excited state at an average 
fidelity of 95.4%. Finally, we showed that the scheme also 
works as a readout with inbuilt classical memory that can 
bring back the measured qubit state at a much later time 
with an average fidelity of 97.3%. The simulations were 
done using a quantum trajectory approach that includes 
fluctuations in the Kerr-resonator and 3d-cavity fields 
and with experimentally realistic parameters. 

The main limitation for the schemes is imperfect tt- 
pulses that suffer from the stochastic perturbations by the 
Kerr-resonator field and finite readout fidelity. Moreover, 
the stabilization scheme suffers from the finite lifetime 
of the qubit, but improving the qubit lifetime and em¬ 
ploying Purcell-filter techniques will improve the fidelities 
significantly. 

Although the intention of this work is to accomplish 
qubit measurement, feedback, and classical memory with¬ 
out extracting information about the qubit state from the 
cryogenic environment, we must have the option to extract 
that information to experimentally verify that the con¬ 
cept works as intended (e.g. for tune-up and debugging). 
In particular, we must verify that after measurement by 
the Kerr resonator, the qubit state and Kerr resonator 
state are properly correlated. This can be established 
with a minimum of added complexity by interrupting the 
protocol in Fig. 4 immediately after the first Ta period. 
At that time, we can measure the state of the qubit by 
applying a pulse into the am port, near resonance with 
the 3d-cavity and of appropriate amplitude and duration 
to constitute a so-called high power readout m- By 
monitoring the reflected pulse phase, even with an in¬ 
efficient microwave ampliher, we achieve a high-fidelity, 
but destructive readout of the qubit state. As the qubit 
cavity system has a much stronger non-linearity than the 
Kerr resonator, and the Kerr resonator is far detuned 
from the qubit cavity at this moment in the protocol, we 
expect that the high power readout can be accomplished 
without altering or being altered by the state of the Kerr 
resonator. Subsequent to the high power qubit readout, 
we can easily determine the state of the Kerr resonator 


by exciting it with a very weak tone added to the ad 
port. The phase shift of this weak tone, monitored again 
through the am port, will depend on the state of the Kerr 
resonator. Because this second measurement detects the 
stable state of the Kerr resonator, it can be performed 
long after the qubit readout when all fields from the high- 
power readout have decayed away. It should likewise be 
independent of the detected qubit state. In this manner 
we can measure the classical correlation between the qubit 
and Kerr resonator state that exist at the time when the 
high-power qubit readout is performed. 

Finally it will be natural to extend the system to include 
more qubits or more resonators. To mitigate the need for 
qubit dissipation for the reinitialization protocol described 
in Sec |IV B] and to suppress residual qubit excitation, a 
second Kerr-resonator can be used to prepare the qubit 
in the ground state, just before the memory state, stored 
in the first Kerr-resonator, is fed back to the qubit. 

As already discussed, the protocols presented in this 
work naturally extend to multiple qubits. Controlling 
many qubits with closed feedback outside of the cryo¬ 
genic environment poses a significant challenge [49| since 
each qubit-readout requires its own signal analyzer and 
a corresponding signal generator to apply the desired 
feedback. For systems with 1000 qubits or more, such 
a setup might be infeasible even with state-of-the-art 
multi-channel analyzers. Therefore, even initialization of 
all qubits may seem an overwhelming task. In contrast, 
the all-cryogenic state-preparation presented in this work 
is inherently parallel, as each qubit has it own memory 
that holds the measured qubit state and autonomously 
controls the input state. The only additional requirement 
to prepare any number of qubit in, e.g., the ground state 
is that enough power is provided and appropriate power 
dividers are installed to drive all qubits sytems. 

A more elaborate extension allows stabilization of e.g. 
a Bell-state of two qubits against phase-flip and bit-flip 
errors using a quantum error correcting code. A possible 
setup for this employs two data qubits in which the Bell- 
state is encoded and two additional measurement qubits to 
detect the errors. Performing two-qubit gates between the 
data and the measurement qubits followed by a readout of 
the two measurement qubits will detect any errors in the 
encoded state [51 [50]. A subsequent feedback on the data 
qubit will stabilize the state. In Fig. [^we have sketched 
how such a scheme might be implemented using a 2d-chip 
with multiple qubits m- Two units, each consisting of 
a Kerr-resonator, a 2d-cavity and a measurement qubit, 
allow us to perform the desired two-qubit gates and the 
additional coupling between the quantum buses and the 
Kerr-resonators will be sufficient to effectively stabilize 
the desired Bell-state in a completely autonomous manner. 
Extending this scheme to a seven qubit code is in principle 
possible and would allow autonomous stabilization of any 
logical qubit state with a limited amount of external 
control lines. 

The Bell-state stabilization represents a minimal im¬ 
plementation of the so-called surface code [501150] . In 
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FIG. 9. (Color online) Schematics to stabilize a Bell-state. The circuit components are the same as in Fig. but instead of 
3d-cavities, we use waveguides on a 2d-chip. The state of interest is encoded into the data qubits, Q1 and Q2, while Ml and M2 
are measurement qubits and they serve the same purpose as the qubit in Fig. To perform the stabilization of Q1 and Q2, we 
apply quantum gates between these and Ml and M2. This can be achieved through the two quantum buses (in yellow). Control 
channels from the Kerr-resonators to each quantum bus allow the feedback to correct errors in the data qubits. 


general, the surface code can be extended to any number 
of qubits, and the outcome of the measurements identifies 
whether an error occurred or not at a given point in the 
code. If an error is found, it is not necessary to correct 
it immediately, but errors can be tracked and corrected 
at the end of the full protocol [50] . All-cryogenics control 
may implement Kerr-resonator memories in a full-scale 
quantum error correction codes and, thus, may become 
vital in future quantum information technologies. 
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APPENDIX: COMPARISON OF MEAN FIELD 
AND FULLY QUANTUM DESCRIPTION 


To compare the semi-classical approach of Eq. (16) 


with the fluctuations of Eq. ( |14| ), we simplify the model, 
such that we only have the Kerr-resonator, that is, we 
set Kf, = 0 in Eq. (16). To describe the Kerr-resonator 
in a fully quantum calculation, we choose to calculate 
quantum trajectories using Monte-Carlo wavefunctions 
(MCWF), such that the unnormalized wavefunction is 
given by the non-unitary Schrodinger equation. 
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FIG. 10. A simulation of the Kerr-resonator with the similar 
parameters as in the main text, but with K(, = 0. The solid 
lines are the semi-classical approach described in Sec. |IIH 
while the dashed lines are the same but using the Monte-Carlo 
wavefunctions of the Appendix. In both types of simulation 
we have averaged over 100 trajectories. The legend indicate 
the value of drive expressed as Od/x/ipi. 


with the Hamiltonian 


_ 

H = HAaO^a + h—a^aa^a + ihy^KdCtdi'i' — «)• (A2) 


Similar to the qubit description in main article, we evolve 
this wavefunction until a the norm becomes smaller than 
a given random number, R, drawn uniformly between 0 
and 1. We then apply the jump operator a to the wave- 
function, we then renormalize and draw a new R. On 
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average this reproduces a full master equation descrip¬ 
tion and it is expected to produce the same behaviour 
as the semi-classical description. Using the unnormal¬ 
ized wavefunction we can calculate any mean value as 
{X) = 1/M with M being the 

total number of trajectories, 

In Fig. we see examples of a MCWF simulation 
and an equivalent semi-classical simulation. We these 
simulation we see that the strong driving and the weak 
driving matches perfectly. We even see that the memory 
effect, which might seem like a classical artefact of the 
semi-classical approach is reproduced in the quantum 
simulation. For the middle driving strength we do, how¬ 


ever, notice a discrepancy between the two methods. We 
ascribe this seemingly different behaviour to a difference 
in the exact bifurcation point. Looking at the individual 
trajectories we see a very similar behaviour between the 
two methods, but for the full quantum approach we find 
10% of the trajectories latching into a high photon num¬ 
ber state, while this number is 22% for the semi-classical 
approach. With these results in hand we can therefore 
conclude that the the behaviour of the semi-classical ap¬ 
proach in the main paper will produce results similar to 
a full quantum implementation, but we expect that the 
required drive strength must be re-calibrated. 
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